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We study the oxygen doping dependence of the equilibrium first-order melting and second-order 
glass transitions of vortices in Bi2Sr2CaCu2 08+5- Doping affects both anisotropy and disorder. 
Anisotropy scaling is shown to collapse the melting lines only where thermal fluctuations are dom- 
inant. Yet, in the region where disorder breaks that scaling, the glass lines are still collapsed. A 
quantitative fit to melting and replica symmetry breaking lines of a 2D Ginzburg-Landau model 
further reveals that disorder amplitude weakens with doping, but to a lesser degree than thermal 
fluctuations, enhancing the relative role of disorder. 



PACS numbers: 74.25.Qt, 74.25.Dw, 74.72.Hs, 64.70. Pf 

Elasticity, thermal energy, disorder, and inter-layer 
coupling are some of the closely competing energy scales 
in the intricate H — T phase diagram of the vortex 
matter in the layered high temperature superconductor 
Bi 2 Sr 2 CaCu 2 8+ a (BSCCO) [J!!!,!!- The low- 
temperature part of the equilibrium phase diagram of 
BSCCO was recently made accessible to experiment by 
vortex shaking 0, mAm- ^ unveiled a first-order (FO) in- 
verse melting line, which continues the thermal melting 
line from high temperatures Q separating low-field or- 
dered phases from high-field amorphous ones. A second- 
order (SO) transition line, at which low-temperature 
glassy phases get thermally depinned, was subsequently 
reported Q. In this letter we study the oxygen doping 
dependence of these transition lines. We show that the 
SO line scales with material anisotropy even where the 
FO line does not, and that effective disorder weakens 
with doping, but gains relative dominance over thermal 
fluctuations. 

We present measurements of optimally doped (OPD), 
slightly over doped (SOD), over doped (OVD) and highly 
over doped (HOD) BSCCO crystals [Tfl, [ll| with critical 
temperatures T c =92, 90, 88.5 and 86 K, respectively, 
corresponding to hole concentrations of 0.171, 0.180, 
0.184 and 0.190 [l§]. Various crystal geometries were 
studied with typical sizes of ^300x300x40 /im 3 . They 
were mounted on 10 x 10 fim 2 Hall sensor arrays, fabri- 
cated in a GaAs/AlGaAs heterostructure. At low tem- 
peratures we utilized a 350 Oc in-plane ac shaking field 
of 10 Hz to relax the pancake vortices towards their equi- 
librium configuration @, @, @|- Conjugating local probes 
with shaking yielded the equilibrium reversible magneti- 
zation of the samples. 

Figure [1] shows the local induction, B(H, T), measured 
in the SOD sample by sweeping the temperature, T, at a 
constant out-of-plane field, H, in presence of an in-plane 
shaking field. A linear term a H T + (3 H was subtracted 
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FIG. 1: (color online) The local induction shows FO (dashed- 
dotted) and SO (dashed) transition lines in the H — T phase 
diagram. The induction in the high-temperature depinned 
liquid and solid phases was artificially flattened by a linear 
subtraction a H T + f3 H (hence the constant color). The low- 
temperature BrG and glassy phases have thus a finite positive 
slope, (a) Values of a H and j3 H used, (b) FO-SO-FO tran- 
sition sequence at a 380 Oe temperature sweep on the back- 
ground of the colorbar. (c) Discontinuous steps in dB/dT, 
manifesting the SO nature of the transition. 



from each temperature sweep (Fig.QJi) to flatten the orig- 
inally increasing (a H > 0) local induction. As a result, 
the FO and SO singular behavior can be readily traced 
throughout the H — T phase diagram. 

The FO melting transition is manifested by a discontin- 
uous step in the magnetization along the dashed-dotted 
line in Fig. [I] It separates the high-field amorphous glass 
and liquid phases from the low-field quasi-long-range- 
ordered Bragg glass (BrG) 0, The melting line be- 
comes nonmonotonic at a certain temperature [H, 0| . Its 
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FIG. 2: (color online) FO and SO transition lines (solid 
and open symbols, respectively) measured with OPD (O), 
SOD(D), OVD (0) and HOD (a) samples. The inset shows 
the doping dependence of the maximum of the FO line and 
its bisection with the SO line. Dotted and dashed lines are 
second-degree polynomial fits. 



inverse-melting part is believed to be an order-disorder 
transition, induced mainly by quenched disorder and not 
by thermal fluctuations 0, 0, [E 0, H, [3 [H, Ell . 

The SO phase transition is manifested by a break in 
the slope of the magnetization at T g , marked in Fig. Q] by 
the dashed line. It separates the fully flattened high- 
temperature magnetization (dB /dT\T>T g — a ~ 0) of 
constant color from a low-temperature region of nonzero 
slope (dB/dT\T<T s — a > 0), whose color varies with 
temperature. The nonanaliticity at T g is demonstrated 
in Fig. [it, which shows ~40 mG/K steps in dB/dT for 
temperature sweeps lying below as well as above the 
FO line. The location of the SO line above melting 
resembles earlier dynamic irreversibility measurements 
17, [HI, 19, 0, I21I and theoretical models of a glass 
H SHU!. Few dynamic mea- 



transition |1|, [! 
surcments found a similar line of thermal depinning be- 
low melting @, [H, 0, 0, ^ , which separates 
the BrG phase below it from a thermally depinned solid, 
whose detailed characteristics are not yet certain. 

We mapped three more samples of various doping lev- 
els that yielded the B—T phase diagram plotted in Fig. [2] 
The melting line shifts as a whole to higher fields with 
doping due to the decrease in anisotropy, which implies 
stronger inter-layer coupling and higher stiffness of the 
pancake vortex stacks [10|, [30|, |3l|, [32 1. As the FO line 



shifts to higher fields, its maximum shifts towards in- 
creasingly higher temperatures (Fig. [2] inset). Identify- 
ing these maxima with the crossover from disorder to 
thermally dominated behavior P, 0, [H, 0, 0| may sug- 
gest that doping enhances disorder. The SO transition 
line also shifts upwards, further signifying that disorder 



is enhanced with doping. However, we will show be- 
low that this conclusion is inaccurate. Moreover, arguing 
naively that the balance between characteristic energies 
constrains the SO line to bisect the FO one at its maxi- 
mum [l[ is also too simplistic. The inset of Fig. [5] clearly 
shows that in over-doped samples the bisection point re- 
sides below the FO maximum temperature (the overlap 
of the two in the OPD sample is apparently accidental). 
Therefore, a naive description of thermal depinning can- 
not account for both the SO line and the FO maximum 
behavior (for instance, the latter is directly affected also 
by elasticity) @,Hli|. 

We now present a quantitative scaling analysis of 
the doping dependence, parameterized by the sample 
anisotropy ratio e 2 = m a b / m c (nii is the electronic effec- 
tive mass in the ith direction). For high-K superconduc- 
tors, such as BSCCO, the Ginzburg-Lanadau (GL) free 
energy functional can be recast into an isotropic (e = 1) 
form by rescaling its parameters [33j |. We focus on one 
such transformation 3J] that rescales space, magnetic 
induction, the penetration depth Ao and the coherence 
length £o- Models for high-temperature melting are usu- 
ally independent of £oj and find Ao to enter with some 
model-dependent power as a proportionalit y fa ctor. For 
definiteness we resort to a specific model [il.l35l| for a FO 
evaporation line with no disorder from vortex solid to 
pancake gas B E (t) cx (e 2 /A§d)(l - t 2 )/t, where t = T/T c 
and d is the inter-plane separation. The scaling transfor- 
mation [2H indeed renders B^t) isotropic. 

Figure [3] shows the rescaled FO and SO lines. The 
high-temperature parts of the FO lines are perfectly col- 
lapsed by dividing their induction axes by a constant 
[301 ] . Bv(t) (dashed-dotted line) fits the collapsed melt- 
ing lines precisely, asserting that the multiplicands in this 
procedure are (e /e) 2 , normalized by the OPD 1 « 500 
[36|. Anisotropy scaling collapses the data down to 
^th ~ 0.58. Below t t h, which appears to be independent 
of anisotropy, the rescaled FO lines disperse again, and 
the fit to Bs(t) breaks. The flattening of the FO line to- 
wards an inverse-melting behavior results from quenched 
disorder, which gains dominance with decreasing tem- 
peratures [l[ . Accordingly, above tth the FO transition is 
purely thermally-induced, and completely unaffected by 
disorder [3]. Just below i t h disorder becomes a relevant, 
though not yet a dominant, energy scale. 

This counterparts the extremely low-temperature be- 
havior, where thermal energy becomes negligible relative 
to pinning, resulting in a flat temperature-independent 
behavior of the FO lines [l|. Indeed the FO lines in 
Fig. tend to flatten towards their ends, below which 
350 Oe - 10 Hz shaking is insufficient for detecting a 
reversible melting step. We thus conjectured a similar 
doping-independent threshold temperature id ~ 0.25, be- 
low which thermal energy becomes irrelevant. We fitted 
the low-temperature order-disorder lines by a leading or- 
der expansion Boo(t > id) ~ -BoD(id) + A(i — td) 2 , which 
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agrees very well with measurement. We can, therefore, 
estimate -Bod (id), at which elasticity is balanced solely 
by the disordering potential. It increases monotonically 
with doping (inset of Fig. [3]), stating that with reducing 
anisotropy elasticity gains dominance also over disorder. 

Yet, the most remarkable outcome of the anisotropy 
scaling shown in Fig.[3]is the simultaneous collapse of the 
SO transition lines (with zero freedom) . The SO lines re- 
side in a temperature region where anisotropy scaling of 
the melting lines fails due to effects of disorder. Still, the 
same scaling transformation somehow succeeds in render- 
ing the glass-transition isotropic, even though this tran- 
sition is believed to be intimately related with the com- 
petition between disorder and thermal fluctuations. 

To gain deeper understanding of the low-temperature 
behavior we fit the measured transition lines to those 



predicted by a recent calculation [37| , which gives access 
to the doping dependence of the model's free parameters. 
It incorporates thermal, disordering and elastic energies 
to yield bisecting FO and SO lines. The pancake vortex 
system is modeled by the 2D GL theory. We interpret 
this single layer model as the outcome of an integration 
of all other layers out of a complete 3D theory, resulting 
in an effective 2D model whose rcnormalized coefficients 
may still depend on anisotropy e of a 3D mass tensor. 

The free energy functional is averaged over gaussian 
disorder in the coefficients of the quadratic and quar- 
tic terms using the replica method. Taking the lowest- 
Landau-level (LLL) approximation yields for the repli- 
cated partition function 

Z*= ( cxpf-^G (* a ) + ^i?|*a| 2 |^ 
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where R = a 2 ^bR/32TT 2 , R is the disorder amplitude in the 
coefficient of the quadratic LLL term, b = B/H C 2, ax = 
(b + t-l)/{2it 2 b 2 t 2 Gi) 1/i is the LLL parameter, Gi is the 
effective Ginzburg number that generally scales thermal 
fluctuations, and k is the GL parameter. Randomness in 
the coefficient of the quartic LLL term, although crucial 
for obtaining the non-analyticity of SO, has negligible 
effect on the locus of the transition lines, and is therefore 
neglected in the present context. 

The FO line is calculated by equating the energies of 
the homogeneous and crystalline states under the influ- 
ence of disorder in gaussian approximation. This ex- 
tends beyond the earlier calculation, which treated dis- 
order perturbatively [38| . The glass line is found from 
the stability analysis of the replica symmetric solution. 
The replica symmetry breaking (RSB) is continuous. The 
corresponding Parisi function describing the hicrarchial 
structure of the glassy state and its detailed derivation 
can be found in Rcf . 13711 . 




FIG. 3: (color online) Data collapse of the over-doped high- 
temperature (t > tth) FO lines (solid symbols) onto the OPD 
line, which also collapses the SO lines (open symbols). The 
fit to E>E(t) (dashed-dotted) sets the anisotropy ratios e/e 
to 1, 1.20, 1.54 and 1.85 for the OPD (O), SOD (□), OVD 
(0) and HOD (a) samples, respectively. The inset shows 
the linear dependence of the characteristic Bod (id) on e/e a 
(dotted line). 



The model's free parameters are R, Gi, T c , and H C 2 = 
T c \ T . However, due to a hidden symmetry, given 

by Gi/H 2 2 and Gi/R 2 , the theory effectively has only 
three independent fitting parameters. We therefore fixed 
H C 2 =100T, for simplicity, consistent with our limited 
range of hole concentrations [l2| ■ Figure [4] presents two 
fitting strategies motivated by different physical behav- 
iors. The first (thick lines) is optimized to fit the col- 
lapsed parts of the phase diagram - SO lines and FO ones 
for t > tth- Interestingly, this set simultaneously provides 
a collapsed FO line also for t < tth- This results from the 
above symmetry of the model in which Gi and e take 
the same role, defining a relation R oc \f(Ti along which 
the calculated phase diagram remains unchanged under 
anisotropy rcscaling. The thick lines in Fig|4] insets show 
the resulting parameters Gi = 3.91 (eo/e — 0.37) 2 and 
R = 1.54(eo/£ — 0.37). The departure of the calculated 
FO line (thick dashed) from the measured ones below 
tth may indicate that theory lacks a symmetry-breaking 
mechanism once anisotropy scaling breaks in experiment. 

An alternative fit (thin lines), faithful to the com- 
plete FO lines, dictates the slightly different set of values 
Gi = 4.77(e /e - 0.38) 2 , R = 0.75(e o /e - 0.2) (insets, 
thin lines). The excellent fit of the FO lines (dashed) 
suggests that the effective 2D model still captures the es- 
sential physics involved. Its validity should break down 
at very low temperatures, where indeed it misses the flat- 
tening of the FO lines as td is approached, and close to T c 
due to critical fluctuations, which may explain the 10% 
overestimate in the fitted bare T c values 391 ]. Substitut- 
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FIG. 4: (color online) Fits of FO (dashed) and SO (solid) lines 
of the LLL 2D GL model with H c2 =100 T, f c =103, 101, 98 
and 95 K for the OPD (O), SOD (□), OVD (0) and HOD 
(A) samples, respectively. The values of Gi (a) and R (b) 
were fixed either by fitting the SO lines and FO ones above 
tth (thick lines) or by fitting the full FO lines (thin lines). 



ing this set of values in the calculation of the RSB lines 
(solid) with zero freedom also produces a good fit that 
improves with doping, but misses the anisotropy rescal- 
ing of the measured SO lines. Nevertheless, the discrep- 
ancy between the two sets of values of the different fitting 
strategies is small, and vanishes in the high doping limit. 

In Fig. the reduced dominance of thermal fluctua- 
tions with doping is clearly captured by the decrease in 
Gi. Its dependence on anisotropy is closer to the 3D 1/e 2 
than to the 2D e-independent behavior. Remarkably, the 
disorder amplitude in Fig. |4]b also decreases with doping. 
This contrasts the result of Ref. (40| . which found that 
oxygen addition increases the defects concentration. The 
apparent contradiction is resolved by noting that the dis- 
order amplitude is affected also by the anisotropy depen- 
dence of the effective pinning potential, which decreases 
with doping as the vortex stacks get stiffer. The over- 
all decrease of R suggests that the latter mechanism is 
dominant, not the rising defect concentration. 

Last, while both Gi and R decrease with doping, Gi 
has a stronger dependence on anisotropy. This clarifies 
the true mechanism that shifts the maximum of the FO 
line and its bisection with the SO line towards higher 
temperatures with doping (Fig. [2|) . Although the magni- 
tude of the disorder amplitude decreases with increased 
doping, disorder still gains relative dominance over ther- 
mal fluctuations, which decrease faster with e. 

In summary, e 2 scaling accounts for the reduced 
anisotropy of BSCCO samples with doping down to ith- 
Below i t h, where disorder becomes relevant, e 2 scaling 
still collapses the SO lines. From the quantitative agree- 
ment with the LLL 2D GL model we deduce that disor- 



der and thermal-fluctuations weaken relative to elasticity 
with increased doping, which shifts the FO line towards 
higher fields. Disorder still gains relative dominance over 
thermal fluctuations, which concurrently shifts the max- 
imum of the FO line and its bisection with the SO line 
towards higher temperatures. 
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